Heat capacity study of BaFe2As2: effects of annealing 
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Heat-capacity, X-ray difTraction, and resistivity measurements on a high-quality BaFe2 As2 sample 
show an evolution of the magneto-structural transition with successive annealing periods. After a 
30-day anneal the resistivity in the (ab) plane decreases by more than an order of magnitude, to 12 
/if2cm, with a residual resistance ratio ~36; the heat-capacity anomaly at the transition sharpens, 
to an overall width of less than K, and shifts from 135.4 to 140.2 K. The heat-capacity anomaly 
in both the as-grown sample and after the 30-day anneal shows a hysteresis of ~0.15 K, and is 
unchanged in a magnetic field /xoH = 14 T. The X-ray and heat-capacity data combined suggest 
that there is a first order jump in the structural order parameter. The entropy of the transition is 
reported. 

PACS numbers: 74.25.Dw,74.25.Ha,74.70.-b,75.30.Fv,75.25.+z 



I. INTRODUCTION 

In the 122 series of Fe pnictide superconductors, super- 
conductivity is produced by doping a parent compound 
AFe2As2 (where A = Ba, Sr, Ca, Eu) with electrons, e.g., 
by substitution of Go^ or Ni^ on the Fe sites or holes, 
e.g., by substitution of K-, Na^ or Cs^ on the A sites. 
It was also found that under pressure the 122 parents 
themselves become superconducting^"— . This led to the 
finding that even substitutions of isovalent elements, for 
example, pi° for As or Ru", Iri^, Pdi^ or Rbi^ for Fe, 
induce superconductivity. 

In this series, the magnetic (spin-density wave) and 
structural (tetragonal to orthorhombic) transitions in 
the parent compound are coincident; there is a com- 
mon magneto-structural transitioni^. The thermody- 
namic nature of that transition, first order or second 
order, remains controversial. In addition, wide ranges 
of the values of the parameters that characterize the 
heat capacity have been reportedi^"— . At least for 
near optimally doped samples, the doping suppresses the 
magneto-structural transition, and the superconductiv- 
ity occurs in the high-temperature tetragonal phase of 
the parent compound. Nevertheless, the heat capacity of 
the low-temperature orthorhombic phase has occasion- 
ally played a role in analyzing the heat capacity of super- 
conducting samples to separate the lattice and electron 
contributions, and thus influenced the interpretation of 
the latter. 

In this paper we report heat-capacity measurements on 
a BaFe2As2 single crystal, as grown and also after vari- 
ous annealing periods. Some preliminary high resolution 
X-ray diffraction measurements are also reported. The 



heat-capacity anomaly at the magneto-structural transi- 
tion evolves in a complicated way with annealing time. 
Although the peak is already relatively sharp in the sam- 
ple as grown, it is substantially sharpened and shifted to 
a higher temperature after a 30-day anneal. In contrast 
with that trend, the heat capacity peak is actually broad- 
ened after shorter annealing times. The sharp transitions 
in the sample as grown and after the 30-day anneal sup- 
port the interpretation of the transition as first order. 
The behavior after intermediate anneals evinces an ex- 
treme sample dependence of its manifestation in the heat 
capacity. 



II. EXPERIMENTAL PROCEDURE 

The measurements were made on a 10.3-mg fragment 
of a large single crystal of BaFe2As2 that had been syn- 
thesized by a modified self-flux methodic. The mea- 
surements were made on the sample as grown, and af- 
ter successive anneals at 700° C under a low pressure 
of Ar gas in a sealed quartz tube, for 1, 4, 8, 14, 
and 30 days. Except for a measurement of the resid- 
ual resistivity, which was made on a sample with a 
well defined geometry, all the measurements reported 
here were made on the same sample. Another piece 
of the same large crystal was used for neutron diffrac- 
tion measurements, which have been reported by Wil- 
son et al}"^ Inductively coupled plasma (ICP) and elec- 
tron microprobe wavelength-dispersive X-ray spectro- 
scopic (WDS) analysis on the as-grown sample con- 
firmed the 1:2:2 stoichiometry. The ICP and WDS mea- 
surements gave, respectively, 1:1. 990±0. 028:1. 995±0. 145, 
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and 1.008±0.05:1.994±0.011:1.998±0.011. X-ray scatter- 
ing attests to the high quahty of the crystal by rocking 
curves having full widths at half maxima (FWHM) of 
about 0.1°. Measurements on the samples were made in 
Quantum Design apparatus: heat-capacity and resistiv- 
ity measurements in a Physical Property Measurement 
System (PPMS), and magnetization in a Magnetic Prop- 
erty Measurement System (MPMS). Synchrotron X-ray 
diffraction data were collected at the 6-ID beam station 
of the Advanced Photon Source, Argonne National Lab- 
oratory. 



III. RESULTS AND DISCUSSION 

Figure 1 shows both the magnetization in a magnetic 
field of 5 T parallel with the (ab) plane (upper window), 
and the zero field normalized resistivity in the (ab) plane 
(lower window) of the as-grown BaFe2As2 crystal. The 
very sharp drop of both the magnetization and the resis- 
tivity at the magneto-structural transition and the small 
low temperature magnetization upturn are signs of the 
high quality of the crystal. 
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FIG. 2: The zero-field C/T of the as-grown BaFe2As2 crys- 
tal. The inset shows the same data near the transition on an 
expanded scale. 
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FIG. 1: Characterization of the as grown BaFe2As2 single 
crystal: magnetization in magnetic field of 5 T // (ab) (upper 
window) and zero field normalized resistivity in (ab) plane 
(lower window). 



The zero-field heat capacity of the as-grown crystal 
is shown in Fig. 2. The magneto-structural transition, 
reflected in the heat capacity by a sharp peak, occurs 
at 135.4 K, where the heat capacity attains a maximum 
value of 320 JK~^mol~^. This value is the highest re- 
ported so far on BaFe2As2 and is another clear indica- 
tion of the high quality of the sample. The inset shows 
the same data near the transition on an expanded scale. 
The "full width at half maximum" of the peak is ^-^0.2 



K, well within the window of possible first-order behav- 
ior deduced from the neutron scattering measurements 
of Wilson et als^ 
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FIG. 3: The low-temperature specific heat of as grown 
BaFe2 As2 crystal. The data, in magnetic fields H // c, are dis- 
placed by progressive increments of 5 mJK~^mol~^ for clarity. 
For each field the curve represents a simultaneous fit to all the 
data. 



The heat capacities of the as-grown sample below 10 
K and in fields up to /LtoH = 14 T parallel with the c 
axis of the crystal are shown in Fig. 3. To the accuracy 
of the measurements, the heat capacity is independent of 
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field, wfiicfi sliows tliat tfie concentration of paramagnetic 
impurities is well below 10"'^ mol/mol. The data in all 
fields were fit simultaneously with the form 



c = 7r- 



Clat-i 



(1) 



where the first term on the right represents the electron 
contribution (Ce) and other terms are the usual approx- 
imation for the lattice contribution (C/at). (The devi- 
ations of the zero-field points below K, which are 
associated with helium adsorbed on the sample and sam- 
ple platform^*, do not make a significant contribution in 
the overall fit.) The derived values of the coefficients are 
given in Table 1, where the uncertainty quoted for each 
parameter is the standard error in the fit. The values of 
B5 and By have a strong dependence on the temperature 
interval and number of terms used in the expression for 
Ciat, and the uncertainties reflects the interdependence 
of those terms. However, the reported values do repre- 
sent the experimental data to within the precision of the 
PPMS measurements, r^2%. Other reported values of 7 
and B3 are included in the table. The experimental data 
are compared with the resulting fit, where the data in 
fields other than zero are shifted in C/T by progressive 
increments of 5 mJK~^mol~^ for clarity. In each case 
the common fitting expression is shown, correspondingly 
shifted. 
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This work 
Ref. 19 
Ref. 14 
Ref. 15 


5.11±0.02 
6.1±0.3 
16(2) 
37 


0.412±0.003 
1.51±0.01 
2.0(5) 
0.60 


8.3±11.0 


6.2±0.9 



TABLE L Coefficients of the low-temperature heat capacity 
(see Eq. 1). The values of the coefficients are in mJ-K-moI 
units. 

In order to investigate the nature of the magneto- 
structural transition, high-resolution X-ray scattering 
measurements of the tetragonal to orthorhombic struc- 
tural transition were performed at the Advanced Photon 
Source's sector 6 beamline. Scattering measurements on 
the 30-day annealed sample were taken with an X-ray 
photon energy of 10.452 keV, in a vertical scattering ge- 
ometry. High resolution reciprocal space scans along the 
[110] direction across the tetragonal (228)t reflection, 
were recorded as the temperature was decreased from 160 
K to a base temperature of 6 K. Close to the structural 
phase transition, measurements were taken at tempera- 
ture steps of 0.1 K. To ensure that the temperature was 
stable during each measurement, a rest period of 4 min- 
utes or greater preceded each scan. 

Figure 4 shows the evolution of the diffraction pro- 
file across the phase transition. Below a temperature of 
139.4 K, two features that correspond to the (408)o and 
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FIG. 4: Reciprocal lattice scans across the (228)t, (408)o, 
and (048) o reflections above and below the structural phase 
transition. The vertical arrows in panels (c) and (d) show 
the (408)o and (048)o reflections. Panels (b) and (c) show a 
continuous splitting of the central (228)t reflection which pre- 
cedes and coexists with the (408) o and (048) o reflections. All 
reflections are fitted with a Lorentzian squared peak profile. 



(048) o orthorhombic reflections appear abruptly. The 
fact that both orthorhombic reflections are visible, indi- 
cates that the sample has orthorhombic domain twinning. 
These reflections (denoted by vertical arrows in Fig. 4 
are independent of the (228)t reflection, and show only 
a gradual evolution with decreasing temperature. They 
are conspicuous at 139.2 K, as shown in Fig. 4(c); they 
had first appeared, with a much smaller amplitude, at 
139.3 K; they were not visible at 139.4 K. This indicates 
that the structural phase transition is first order in na- 
ture. Surprisingly however, this phase transition is pre- 
ceded, at 140.1 K, by a splitting of the (228)t reflection, 
which progresses continuously from zero with decreasing 
temperature. At 139.6 K the splitting still appears as a 
broadening of the (228)t reflection, which can be decom- 
posed into two components, as shown in Fig. 4(b). At 
139.2 K it has increased to the point that two reflections 
are clearly visible. At lower temperatures the amplitude 
of this reflection decreases rapidly and the splitting con- 
tinues to increase. This reflection is barely visible at 
139.0 K, and not at all at 138.7 K, as shown in Fig. 4(d). 
It seems probable that the two components approach the 
orthorhombic reflections, but a quantitative treatment 
of the kind shown in Figs. 4(b) and 4(c) is not practi- 
cal. In any case, it is noteworthy that both components 
of this central reflection disappear as the orthorhombic 
reflection increase, as expected from a first order tran- 
sition between two ordered phases. The data in Fig. 4 
were all obtained on cooling but the major features are 
reproduced on warming. This is demonstrated in terms 
of the structural order parameter in Fig. 5, where both 
cooling and warming data were included. The very small 
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discontinuities in the order parameters at 141 K (in the 
region between Figs. 4a and 4b) are not significant in 
relation to the uncertainty in the decomposition of the 
central reflection into two components. 

Therefore this result indicates that there are in fact two 
parts to this phase transition. First, below 140 K, the 
tetragonal structure becomes orthorhombic with a small 
distortion. The x-ray scattering data shows a splitting of 
the tetragonal reflection into the two orthorhombic reflec- 
tions, that starts continuous with temperature. Second, 
at the slightly lower temperature of 139.4 K, a transition 
between two orthorhombic structures is observed. This 
change between the two orthorhombic structures hap- 
pens abruptly with temperature. Thus, these results are 
consistent with the thermodynamic measurements that 
shown a first order transition. To emphasize this point we 
present the structural order parameter, (5=(a-b)/(a-l-b), 
as a function of temperature in Fig. 5. The lattice con- 
stants a and b were determined through fitting the posi- 
tions of the reflections with a Lorentzian squared profile 
in a non-linear least squares analysis. Figure 5 shows 
that the two sets of orthorhombic reflections co-exist for 
a temperature range of 0.4 K. It should also be noted 
at this point that the ordering temperature of the first- 
order transition is slightly lower than the temperature 
indicated by the heat capacity measurements. Possible 
explanations include slight differences in the calibrations 
of the respective thermometers and, for the X-ray ex- 
periment, a subtle heating of the sample by the X-ray 
beam which is not recorded by the temperature sensors. 
It should also be noted that the stability of the temper- 
ature during each scan and the reproducibility of these 
results upon cooling and warming, indicate that a kinetic 
origin of these effects is unlikely. 
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FIG. 5: The structural order parameter, (5=(a-b)/(a-|-b), as a 
function of temperature. Filled symbols are for cooling data 
and open symbols are for the warming data. 



The unusual behavior of this structural transition may 
provide an explanation as to why previous diffraction 



measurements have characterized it as second order. 
Clearly, with lower angular and temperature resolutions, 
one could mistake the initial splitting of the tetragonal 
reflection and the discontinuous appearance of the or- 
thorhombic reflections as being the same thing. In fact 
Wilson et alM- acknowledge this point by stating that 
the phase transition could be first order within a temper- 
ature interval of 0.5 K. (Their measurements of the struc- 
tural order parameter, by neutron scattering, were made 
on another piece of the crystal from which the as-grown 
sample was obtained. They reported a temperature de- 
pendence of 5 similar to that shown in Fig. 5.) Whether 
the effects described in this paper are a universal feature 
of this phase transition, or just specific to this annealed 
sample, remains to be determined. However, these re- 
sults do make a compelling case for further investigations, 
especially their relationship to the concurrent magnetic 
phase transition. Specifically, simultaneous high resolu- 
tion diffraction measurements of both the magnetic and 
structural order parameters will be necessary to clarify 
the precise nature of the phase transition(s). 
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FIG. 6: (a) Nomalized resistivity of BaFe2As2 single crystal, 
as grown and after 1, 4, 14, and 30 days of annealing, (b) The 
same data as in (a) near the transition. The data were shifted 
in p(T)/p(300 K) to facilitate comparison of the transition 
temperature after annealing, (c) Residual resistivity ratio of 
in-plane resistivity, RRR=p(300 K)/p(2 K), versus annealed 
time periods. 



Figure 6a shows the temperature dependence of the 
resistivity in the (ab) plane normalized to the resistiv- 
ity at 300 K for the as-grown sample and after anneal- 
ing periods of 1, 4, 8, 14, and 30 days. The transition 
temperature increases continuously with annealing of up 
to 4 days, but remains constant after longer annealing 
periods (Fig. 6b). The resistivity has a weak linear de- 
pendence on temperature above the transition but de- 
creases sharply below the transition. This is in disagree- 
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ment with the log{l/T) behaviour below 12 K reported 
by Wang et alJ^ In Fig. 6c the in-plane residual re- 
sistivity ratio, RRR=p(300 K)/p{2 K), is plotted versus 
annealing time. While the RRR of the as-grown sample is 
slightly above 5, one of the highest value reported for the 
system, it is continuously increased by annealing, attain- 
ing a value of 36 after 30 days. The residual resistivity, 
p{2 K), of the sample after the 30-day anneal is about 
12 /iJlcm, one order of magnitude smaller than the low- 
est value reported so far on the system. This decrease 
of the residual resistivity can be caused by the healing 
of the imperfections and dislocations that are centers of 
scattering through the annealing process. 
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FIG. 7: C/T near the transitions of the BaFe2As2 sample, as 
grown and after annealings. 



The heat capacity near the transition is shown in Fig. 
7. After 4 or 8 days of annealing the single sharp peak 
of the as-grown sample evolves into a broad feature at 
a somewhat higher temperature with an overall width 
of ~1.5 K. With a longer annealing time, 14 days, the 
anomaly moves to still higher temperature, with approx- 
imately the same width. After 30 days of annealing the 
broad feature evolves back into a single sharp peak at 
140.2 K, with the heat capacity reaching a maximum of 
600 JK^^mol^^. The entropy change at the transition of 
the sample, both as grown and after the 30-day anneal is 
-0.84 JK-^mol-i. 

The evolution of the heat-capacity anomaly at the 
transition with annealing time is unusual. Annealing a 
sample, at least under suitable conditions, is usually ex- 
pected to produce a steady improvement in sample qual- 
ity. The resistivity results are consistent with that ex- 
pectation. The behavior of the heat-capacity anomaly 
- broadening at intermediate annealing times, followed 
by sharpening at longer times - has no obvious expla- 
nation. It does, however emphasize the sensitivity of 
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FIG. 8: C/T near the transitions of the BaFe2As2 sample after 
the 30-day anneal in zero field, heating (■) and cooling (♦); in 
14 T, heating (•). The inset shows the relaxation curves of the 
sample temperature versus time of a measurement spanning 
the transition. The arrows point to regions of the curve that 
have a reduced dT/dt, showing the presence of the latent heat, 
a signature of the first order transition. 



the transition to details of the sample preparation. It 
is also consistent with measurements on other samples 
that have shown a variety of heat-capacity anomalies as- 
sociated with the transition. The reported widths of the 
specific-heat anomalies, at their "base" varies from 1 to 5 
K. A variety of peak heights, in C/T measured from the 
base, and transition temperatures have been reported: 
0.3 JK-^mol-i at 138 ¥M, 0.9 JK-^mol"! at 136 K^^, 
0.1.5 JK-^niol-i at 132 K^^, 0.4 JR-^mol-i at 133 
K3i, and 0.5 JK^^mol"^ at 133 K^^. In several cases 
the anomalies are comparable to those reported here for 
intermediate annealing times. In one case^S. a double 
peak was reported, i.e., similar in that respect to those 
reported here for intermediate annealing times. 

Figure 8 shows C/T in the vicinity of the transition 
after 30 days of annealing. These data (and the corre- 
sponding data in Fig. 2) were obtained by an analysis 
of the heating and cooling segments of individual heat- 
ing cycles, with a program that was provided with the 
PPMS for use in the vicinity of sharp features in the 
heat capacity. The zero-field data derived from the heat- 
ing and cooling segments show a hysteresis of 0.15 K. The 
existence of hysteresis is consistent with the conclusion, 
based primarily on other results, that the transition is 
first order. (This subtle hysteresis of 0.15 K would not 
have been detected in the neutron study of Wilson et 
al—-) The inset of Fig. 8 shows the heating and relax- 
ation curves for a heat capacity point that spanned the 
transition. The arrows point to regions of reduced dT/dt, 
which show the presence of a latent heat, another indica- 
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tion of a first-order transition- The in-field data show 
that the transition is unaffected by a magnetic field iiqH 
= f4 T. This is a surprising result given the combined 
magnetic and structural nature of the phase transition. 
The heat capacity of the as-grown sample in the vicinity 
of its transition showed the same evidence of hysteresis, 
a latent heat, and the absence of an effect of a magnetic 
field. 
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FIG. 9: Zero-field, low-temperature C/T versus of the 
as-grown and 30-day-annealed BaFe2As2 samples. The solid 
curve represents the collective fit to the data for the as-grown 
sample in all fields. 



The zero-field low-temperature heat capacity after the 
30-day anneal is compared with that of the as-grown sam- 
ple in Fig. 9. The measurements on the annealed sample 
were made on '^1/2 of the sample, which accounts for 



the scatter in the data. To within the precision of the 
measurements, the parameters that characterize the low 
temperature heat capacity (see Table I) are unaffected 
by the anneal. 



IV. SUMMARY 

We made heat-capacity. X-ray-diffraction, and resis- 
tivity measurements on a high quality BaFe2As2 sample 
as grown and after annealing for 1, 4, 8, 14 and 30 days 
at 700° C. We observed an evolution of the magneto- 
structural transition with successive annealing periods. 
After 30 days of annealing the resistivity in the (ab) plane 
decreased to 12 /Ltflcm, with a residual resistance ratio of 
^36; the heat-capacity anomaly at the transition sharp- 
ened, to an overall width of less than 1 K, and shifted 
from 135.4 to 140.2 K. The heat-capacity anomaly in 
both the as-grown sample and after the 30-day anneal 
shows a hysteresis of 0.15 K, and is unchanged in a mag- 
netic field ^qH = 14 T. The X-ray and heat-capacity 
data combined suggest that there is a first order jump in 
the structural order parameter. This jump is also appar- 
ent in regions of reduced dT/dt in the heating and re- 
laxation curves taken in the heat-capacity measurements 
near the transition. The conduction-electron density of 
states and the low-temperature lattice heat capacity of 
the orthorhombic phase are unaffected by the annealing. 
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